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ABSTRACT 


A new remote sensing technique was developed for pre- 
dicting the expected mean laSer radar return froma rough 
water surface. This technique involved measuring the stan- 
dard deviations of the upwind and crosswind profiles of the 
elliptical glitter patterns occurring for illumination of 
the water surface with a point source near the laser radar 
system. A pencil beam laser radar from a companion project 
Simultaneously meaSured the reflected signals from the water 
Surface. The glitter pattern images were recorded with a 
Video camera and recorder. The images for each run were 
later digitized along their major and minor elliptic axes 
and averaged over 256 images to produce smooth intensity 
curves from which the standard deviations were measured. 
The radar return fluctuated over a large range because of 
the rapid variation of individual water surface facets, and 
SO waS recorded and time-averaged over the same interval as 
the video images. Data sufficient for empirical prediction 
of expected mean laser return Signal were obtained. This is 
necessary to permit evaluation of the performance of a given 
laser radar design. The data obtained also approximated the 


predictions of a new model proposed in this work. 
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A. BACKGROUND 

The situation at the outset of this project was that a 
working model of an inexpensive, light, and compact laser 
radar altimeter, for use over the ocean surface, had been 
developed in this laboratory in a companion project 
[Ref. 1]. However, it had not been possible to relate its 
performance in the laboratory to real field performance with 
much certainty. The device determined the height above the 
water by measuring the time delay between the emission of a 
short laser pulse and the arrival of a return echo from the 
water surface. The limitations of size, weight, electrical 
power, and cost dictated a low power diode laser system, 
Additionally, there was a requirement of a finite divergence 
angle of the laser beam and associated detector optics in 
order to assure triggering of the time-delay measuring sys- 
tem in the event of a glassy calm and some angular wobble in 
the attitude of the laser system. These collective require- 
ments meant that the primary design problem was that of 
achieving sufficient return signal from the water to assure 
triggering of the time delay circuits. 

Evaluating the performance of the existing working model 
of a laser radar altimeter system over the water had the 
problem that the reflectance of a water surface is dependent 
on the roughness state of the water. The only data in the 


literature known to this author, which would permit an 


estimate of the laser reflectance of rough water, are those 
of Petri [{Ref. 2]. These data are rather meager, consis- 
ting of 16 separate measurements, and requiring a know- 
ledge of the wind speed. The latter had been measured at a 
height of 60 feet above the watCere The wind speed at the 
water surface was thus quite uncertain. Other data by Cox 
and Munk [Ref. 3] related the glitter pattern profile, for 
reflection of the sun, to wind speed over the ocean. Their 
optical data were taken from an aircraft at 2000 feet alti- 
tude, with windspeed measured on a ship at two heights, 9 
feet and 41 feet above the water. It was not clear which 
of these two heights was used for quoting the wind speed. 
Their data also were not directed at evaluating the magni- 


tude of the reflectance, 


B, GENERAL OBJECTIVE 

Preliminary field experiments, carried out as part of 
the work reported here, indicated that the reflectance of 
rough water could vary rapidly under changing wind condi- 
tions, and that knowledge of the wind speed at a given 
instant did not serve as a good indicator of the optical 
properties of the water surface at that instant. Addi- 
tionally, if wind speed were to be measured, it probably 
Should be measured very close to the water surface, This 
would be difficult, expensive, and in many cases impracti- 
cal. The field work reported here for this project was 


carried out from high bridges, but it was intended that the 


work would be extended later to aircraft. In that case the 
necessity for associated ship measurements of wind speed 
would be difficult logistically. Consequently, the objec- 
tive of this work became to develop and demonstrate the 
feasibility of a remote sensing. technique for evaluating 
the reflectance of a rough water surface at the instant of 
a laser radar altimeter test. 

The papers of Cox and Munk, and Petri, taken together, 
led to a conclusion that it might be possible to predict an 
expected laser radar return on the basis of a measurement 
of the width of glitter patterns. It also seemed likely 
that glitter patterns measured with a nearly point source 
of light at night would be simpler to obtain and would 
yield more directly the desired result than would glitter 
patterns obtained with the necessarily oblique illumination 
from the sun. Consequently, the primary objective became 
to test the feasibility of prediction of laser return from 
the profile widths of glitter patterns obtained with a 
point light source. The early tests indicated that laser 
return was in fact a function of glitter pattern width. A 
tentative theoretical model was developed, which, although 
it has not as yet succeeded in predicting the observed 
absolute magnitude of laser return, does lead to a func- 
tional relationship which is useful. With the present 
semiempirical relationship, laser return can be related 
quantitatively to Lambertian return. Although more experi- 


mental data is needed to refine the results, and further 


10 


theoretical work is needed to clear up the discrepancy in 
the absolute values predicted by the model, laser return 
signals measured in the laboratory with Lambertian targets 


can now be directly related to field-expected values. 


C. ASSOCIATED PROBLEMS 
1. Laser Signal Fluctuation 

There were several additional problems in carrying 
out the field experiments at the outset of this project. 
The laser radar altimeter model to be tested used a fairly 
narrow laser beam of 1.2 by 2.4 degrees divergence, with a 
3 degree circular cone of acceptance of the receiver op- 
tics. This divergence had been chosen in the design pro- 
cess of the companion project as an optimization of the 
transmitter-receiver optics under the confines of required 
return signal magnitude and availability of inexpensive 
commercial components. This narrow beam was incident on 
the center of the "glitter" pattern, about which much more 
will be said later. The beam will be treated here as a 
pencil beam, essentially evaluating the peak value of laser 
return at the center of the glitter pattern. A consequence 
of the narrow beam geometry was that reflection of the 
laser beam from the water surface occurred at a relatively 
small and finite number of glint facets. As a result, the 
laser return signal fluctuated rapidly with time over a 
wide range of magnitude. This required recording the 


signal for later data processing to obtain the mean and 


1] 


standard deviation of the variation. In order to accom- 
plish the recording it was also necessary to provide a 
pulse stretcher because the laser return pulses were other- 
wise too short to handle with most analog FM data recorders, 

Although the rapidly fluctuating magnitude of the 
laser return posed a measurement problem, it is helpful in 
the practical application of a laser radar altimeter. Mea- 
surement of range is accomplished by determining the time 
delay between the primary pulse and the return pulse. With 
fairly rapid pulse repetition rate, an occasional high 
pulse will serve to make the range determination practical, 
where the mean signal may not be above threshold for satis- 
factory operation. This feature tends to favor a small 
divergence beam. The statistics of this were not analyzed 
here but are left for later work. The data tape recordings 
from the field experiments are now available for such 
analysis. 

2. Calibration by Means of Specular Reflection 

Prior to the work reported here, optimization of 
design and laboratory evaluation of the laser radar altime- 
ter had been carried out primarily using specular reflec- 
tion from the outside of a plate glass window to obtain a 
return signal, after careful measurement of the inevitable 
slight surface curvature of the window, This permitted 
development of the time-difference circuitry and avoided 
the difficult problem of obtaining a water surface at 


sufficient range to be meaningful. The constant magnitude 
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of the return signals in this situation was also helpful in 
the development work, In evaluating the magnitude of the 
signal, the difference of reflectance of the two glass-air 
interfaces from that of a single water surface was easily 
accounted for. A rough evaluation of system performance 
was obtained in this manner, but the return signal, if 
maximized by angular adjustment to obtain the maximum sig- 
nal, evaluated only the peak value of the central maximum 
of the laser beam, as will be discussed in Appendix B. To 
obtain a measure of the expected return signal from a rough 
water surface or a diffuse reflecting surface from the 
specular reflection data required carrying out a two-dimen- 
sional integration of the laser beam flux over the laser 
beam cross-section profile. The beam tended to have consi- 
derable detailed "banding" structure in its profile, making 
such evaluation rather doubtful. 
3. Calibration by Means of Lambertian Reflection 

The problems attendant on use of specular reflec- 
tion from a plate glass window to obtain a return signal 
are alleviated by use of a diffuse, "Lambertian" reflector 
such as a white sheet of paper. For a finite area, such a 
surface reflects a power proportional to the cosine of the 
angle of incidence, Such surfaces were used in the labora- 
tory. However, the Lambertian surface must be large enough 
to inelude the entire transmitted laser spot. Hence at 
ranges sufficient to simulate real field conditions, the 


area required was large enough to make obtaining a good 


i 


uniform reflectance white surface difficult. The signals 
with such a target would be constant in magnitude, except 
for a small amount of atmospheric scintillation. At close 
range, the signals were so large as to make circuit devel- 
opment for small signal conditions difficult. Calibration 
of the return signal magnitude can be done with a Lamber- 
tian surface at close range, by careful attention to quan- 
titative attenuation of the large return signal that would 
otherwise saturate the receiver circuits. This sort of 
calibration was done in the field tests to be reported 
here. Finally there seemed to be no fully satisfactory way 
to simulate all the circumstances to be encountered by a 
real life laser radar altimeter, other than by field tests 
over wind-roughened water. These were then initiated and 
the accompanying remote sensing techniques reported here 


Carried out, 


D. GENERAL METHOD 

This work reports a new remote sensing system for esti- 
mating the expected laser radar return signal through an- 
alysis of glitter patterns. This is a remote sensing 
method permitting determination of the expected reflectance 
of a rough water surface using equipment colocated with the 
laser radar altimeter. The equipment for this technique 
also offers the advantage of being inexpensive and phy- 


sically small for field use, 
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The time-averaged statistical distribution of light 
intensity was measured for images of glitter patterns pro- 
duced by reflection of light projected onto the water from 
a small area incandescent source located near the laser 
altimeter. The images of the glitter patterns were re- 
corded with a video camera, which avoided the nonlineari- 
ties of photographic recording and permitted digitization 
without intermediate steps. The magnitude of the returned 
laser signal was recorded at the same time that the video 
images were obtained, Because this signal was rapidly 
varying, it was recorded for statistical processing later, 
The resulting data related the mean laser return to the 
glitter pattern profile parameters. The results are most 
succintly expressed as a ratio of reflectance of the rough 
water to the reflectance of a (virtual) Lambertian reflec- 
tor which reflects the same fraction of the total light as 
a water surface (0.0204) at normal incidence. Through use 
of data in other references, the return was also related to 
an inferred wind velocity for comparison. An analytical 
model was developed which seems to predict the correct 
functional behavior, although the predicted absolute magni- 
tude of reflectance is at present out of line. Although 
these experiments should be repeated for refinement and 
completion, these results allow evaluation of a practical 
laser radar system by reflection from a Lambertian reflec- 
tor within the laboratory, providing a more immediate and 


readily controlled environment for equipment development. 


lo 


I] THEORETICAL CONSIDERATIONS. 


A. GENERAL 

When light from a small diverging source is incident on 
a water surface roughened by wind driven waves, it is re- 
flected from the water surface in a glitter pattern, as 
seen by an imager located near the source, The glitter 
pattern is approximately elliptical in shape, with the 
major axis of the ellipse in the upwind direction, and the 
minor axis in the crosswind direction [Ref. 4]. The pat- 
tern is produced by specular reflection at a large number 
of rapidly changing facets located at the points where the 
surface is perpendicular to the incident light. A repre- 
sentative glitter pattern is shown in Figure 1. 

Laser reflectance from calm waters can be readily cal- 
culated analytically [Ref. 5]. This is the case of specu- 
lar reflection whereby a mirror image of the source is 
created equidistant below the water surface, Reflected 
energy appears to originate from this source, 

At the other extreme from specular reflection is dif- 
fuse reflection. In this case a surface scatters incident 
energy over all angles in a hemisphere, resulting in a 
distribution called Lambertian, in which the reflected 


energy varies as the cosine of the angle from the perpen- 


dicular, 
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beetire ¢. Glitter pattern produced by wind speed 
of approximately 12 miles per hour. 


tg 


Between these two extremes, an optically rough water 
surface exhibits properties of both types of reflection, 
The reflected intensity distribution varies according to 
the distribution of wave slopes scattering the incident 


beam, [Ref. 6]. 
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B. GLITTER PATTERNS 

The technique developed in this work was a method of 
evaluating the expected laser radar return from a water 
surface when it is roughened by wind. The technique is a 
remote sensing technique in that all measurements are made 
from the general location of the laser radar altimeter. 
Thus it does not involve knowing the windspeed at the water 
surface, nor any other meteorologie or oceanographic par- 
ameters. The measurements involve only the optical proper- 
ties of the water surface, This is done by determining the 
size of the glitter pattern from a video camera recording 
of the reflection of an approximately point source of light 
near the laser radar. The laser radar used a pencil beam, 
so that it responded to the central maximum reflection 
point of the glitter pattern profile. The reflection from 
this fluctuated through a large range, as a function of 
time, so that a time-average of a recorded signal was 
required. 

Cox and Munk have developed a method of measuring water 
surface roughness from photographs of sun glitter [Ref. 3]. 
The wave slope distribution was deduced from observations 
of sun glitter patterns, with the assumption that the wave 
slope distribution and the reflected light intensity dis- 
tribution were the same. They found that a Gram-Charlier 
distribution was the best fit to the experimental data. 


This distribution function is 
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upwind/downwind direction, and 

Y = @O/o, , with @ the angle from the vertical 
in the ecrosswind direction, 
and o 


The quantities o @ are the upwind and crosswind 


u 
"standard deviations" for this distribution, For a Gaus- 
sian distribution the standard deviation is the root mean 
square deviation from the mean. The first two terms of the 
Gram-Charlier distribution are the same as for a Gaussian 
distribution, The Gram-Charlier sigma values are the val- 
ues of sigma in the second term, As will be discussed 
later, the peakedness of the Gram-Charlier distribution 
causes these sigma values to be located lower on the curve 
than for a Gaussian distribution, The terms in the brac- 
kets in the Gram-Charlier distribution are a form of Her- 
mite polynomial. These terms express the peakedness and 
Skew of the function, 

The Gram-Charlier distribution, above, is a two di- 


mensional distribution, It is displayed as a function of 
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6 and @ in Figure 2. The dotted curves are the shape of a 
Gaussian with the same standard deviation values, The 
crosswind distribution is symmetric about the vertical but 
is more peaked than Gaussian. The upwind distribution is 
also more peaked than Gaussian, and has a skew in the 
upwind direction. 

The Gram-Charlier distribution is a normalized func- 
tion, i.e., with the complete Hermite polynomial series 
included, the integral of the function from minus infinity 
to plus infinity in X, and from minus infinity to plus 
infinity in Y, gives unity. At first sight it would seem 
that integration from minus infinity to plus infinity, over 
the angular variables used, is not strictly correct, as the 
angles can only go to x/2. The function as written above 
is an approximation, as it has only the first few terms of 
the Hermite polynomial ineluded, However, Cox and Munk 
stated that this approximation was valid for all their data 
out to X = Y = 2.5. Contributions to the integral for 
values of X and Y beyond 2.5 are negligible, so that writ- 
ing the integral as if it extended to infinity should be of 
no concern. The largest sigma values encountered in the 
data reported here were 7.5 degrees. For X = 2.5 this 
means the expression is taken to be valid to 19 degrees, 
The integrals written later in this thesis will express the 
limits as minus infinity to plus infinity, recognizing that 


the contributions of high X and Y are negligible. 
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Figure 2. Wave slope distribution profiles deduced by 
Cox and Munk [Ref. 3]. Upper curves are taken along the 
Minor (Crosswind) axis of the elliptical pattern. Lower 
curves are taken along the major (up/downwind) axis. Solid 
curves are distributions deduced for a wind speed of 10 
meters per second, Dashed curves refer to Gaussian distri- 
butions of the same standard deviation. 
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The extent to which the Gram-Charlier distribution 
exceeds the Gaussian due to the peakedness at the center 
can be evaluated by inserting the values of the coef- 
ficients, as given by Cox and Munk, and letting X = 0 and Y 
= Q. For the crosswind direction, the curve has no skew 
and the peakedness is independent of wind velocity. The 
coefficients for crosswind peakedness are 


C40 = 0.40 + 0.23 , 


i! 
& 
e 
=) 
wt) 


C4 £0.41", (2) 
and Coo = 0.12 + 0.06. 
This gives a term multiplying the value of the peak as 
given by the first two terms of equation (1), i.e. the 
height of the Gaussian central peak is multiplied by a 
factor f = 1.109 + 0.061 . (3) 
This factor will be used as 1.11 later. 
For the upwind direction the skewness constants are 
Co, = 0.01 - 0.0086W + 0.03 , 
and Cg3 = 0.04 - 0.033W + 0.12 , (4) 
where W is the wind speed in meters per second. Although 
these constants depend on wind speed, they do not affect 
the area under the curve, and hence the normalization, as X 
and Y appear only in odd powers in the Gram-Charlier Hermi- 
tian terms. Thus whatever area they add for positive 
values of X or Y, they subtract for negative X or Y. 
Although the skewness terms do not affect the peaked- 
ness for the crosswind direction, they do make a small 


contribution to the peakedness in the up/downwind 
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direction, but the contribution is small. For example, for 
a wind speed of 10 meters per second, as illustrated in the 
lower curve of Figure 2, the additional peakedness is 
approximately two percent. The position of the peak is 
also shifted slightly upwind, However, its position is 
still within the "pencil" beam of 1.4 x 2.8 degrees used in 
the laser radar altimeter here. The values for the predic- 
ted laser return signal presented later in this thesis have 
been adjusted for the skewness contribution to the peaked- 
ness, even though this adjustment is in most cases the 
order of one percent - really negligible in view of the 6 
percent uncertainty quoted above from Cox and Munk on the 
value of the peakedness factor, 

To return to the meaning of the "standard deviations" 
for the Gram-Charlier distribution, the standard deviation 


values, 9, and oa are the standard deviations of the 


Cc u? 


Gaussian term of the Gram-Charlier form. In measuring 
these quantities on an actual measured distribution, as 
will be described later in this work, it is necessary to 
measure the distribution curve width at a point lower by a 
fraction 1.11 (or 1.12, depending on the skew contribution 
to the peakedness for that particular curve) than the usual 
0,659 times the peak value - the case for the usual Gauss- 
ian distribution. This has been done in evaluating the 
data later in this thesis. 

A further fine tuning should also be noted. The peak- 


edness terms in the Gram-Charlier distribution do produce 
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upward, or downward, displacements of the curve throughout 
the curve. This can be seen in the curves presented in 
Figure 2. Because of this displacement, it might seem that 
determination of sigma values by meaSuring the width of the 
peak at a height equa] to 0.659/1.11 times the peak value, 
would be in error. (The factor 1.11 comes from equation 
(3) above.) To evaluate this effect, the coefficients in 
equation (2) above produce vertical displacements of the 
Crosswind curve, i.e. for Y = 0 , at X = +1 (at the o, 
point), given by 
ome Ge Cle = ote SY 
+ (1/24) Coy (3) 
=e eae see 0192 
+ .0280"+ .0513 
= eeer004Set .0547. 

This represents a downward shift of 0.4 percent, with an 
uncertainty of 5.5 percent. The 0.4 percent correction 
could be made by drawing a parallel line to the curve, 
shifted downward by that fraction, but it seems not worth 
doing in view of the 5.5 percent uncertainty. For the 
up/downwind curve the corrections are equally insignificant. 
Consequently such corrections have not been made during 
data reduction. 

PiemmOct wor miticant Quantity, Lor this work, is the 
ratio of the peak of the distribution to the area under the 
distribution, as a function of the upwind and crosswind 


Standard deviations. As mentioned before, the quantity, G, 


25 


in equation (1) above, is a statistical distribution. It 
is normalized so that the area under the curve is equal to 
unity, 1,e. 
ad a ~(X%4+y2)/2 
Area = l dea (i etea elie ([H]) dk dYsaaeoy 


where the Hermitian peakedness and skewness terms have. been 


abbreviated by the bracket [H]. 


The power per unit solid angle, I, reflected from a 
small spot on the water surface that is illuminated by the 
laser radar beam will have its reflected light spread into 
a cone where the intensity distribution varies according to 
the Gram-Charlier distribution. The magnitude of the solid 
angle of the cone will be taken to be that subtended by the 
glitter pattern, as observed at the detector. This assump- 
tion is discussed in Appendix D. 

Letting I, be the power per unit solid angle reflected 
perpendicular to the water surface and hence back to the 
detector of the laser radar, the power per unit solid angle 
in any direction will be given by 

I= 1, {2xo,¢,/(H,1]} G , (6) 
where [H,J] = 1.11 for peakedness at X = 0 and Y=0. 

As a verification of the statement in equation (6), 


note that at KX = 0 and Y = Q 


G=G, = {1/(2ne,6,)} [Hy] . 


Hence at X = 0 and Y = Q 


1 = 1,4 (2x .0,)/ 1H, 1Ve( 8) (2x ce) see 
and I =I, ; as tt should be. 


Now the total reflected power is given by 
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P = bf le XeeGiy oo. 


“co =co 


Inserting I from equation (6) gives 


P = If 1, (2n0,¢,/[H]) G dX dY . 


-Cco -& 


Then from equation (5) 
ee A ADS Os ee (7) 
Now, the total reflected power, P, is the power radiated by 


the lasar radar in the pencil beam, P multiplied by the 


0? 
Fresnel reflectance of water at normal incidence, r , or 
Peoegeeeset lt. {2nc.o/[H,)) (8) 
Inverting gives 
Penh | (0/ 2nc, ceeevl |e. 3) 
The flux at the detector, Fog» the power returned per unit 
area of the detector receiving optics, is given by I, times 
the solid angle subtended by a unit area at the detector 
Optics, or 
x 1/R? 


E I 


og O 


Pie. (WM Cie a [ists Ie ae (10) 


C. FRESNEL REFLECTANCE 

Use of the water reflectance for normal incidence in 
the preceding section perhaps needs some discussion. Dur- 
ing the measurement of a glitter pattern the reflection at 
each glint is at normal incidence, Similarly, for the 
light seattered back to the laser detector from a pencil 
laser radar beam, the angle of incidence is again zero. On 
the other hand, for the light that does not return to the 


receiver because it is reflected out into an elliptic cone, 


2] 


the angle of incidence at the reflecting facets is not 
zero. However, the contribution of scattered light for 
facets with angles of incidence greater than 2 sigma is 
negligible. Since the maximum sigma encountered in this 
work was 7.5 degrees, the variation of reflectance for 
angles of incidence of this magnitude needs to be examined. 

The reflectance of a dielectric surface of index of 
refraction, n, as given by the Fresnel equations, for light 
with its electric vector perpendicular to the plane of 


incidence is 
cos © - (n2 - sin@ @) 1/2 4 


ate = (11) 








cos @ + (ny - sin? Q)1/2 
For light with its electric vector parallel to the plane of 


incidence, the reflectance is 








-n2cos @ + (n2 - sin? @)1/2 12 ; (12) 
. =) SS eee é 12 
P n2cos © + (ny - sin? @)1/2 
At normal incidence, © = 0 andr, = rh =f. 


For unpolarized light at other angles, the total reflec- 
Cane mi Sint ae ry )/2 . The reflectance, r, at some repre- 
sentative angles, for water of index n = 1.333, at 16C, is 


r r r Factor 
Normal incidence De 02037 02057 .02037 1.000 
oe .02133 .01944 .02038 1.001 
20° .02449 .01662 .02055 1.009 
30° .03093 .01194 .02144 1.053 
40° .04316 .00585 .02450 1.203 
45° .05299 .00281 .02790 1.370 
Brewster's angle 53.1° .07817 .00000 .07817 3.838 
Grazing incidence 90° 1.0 


As can be seen from the above, the total reflectance for 
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unpolarized light varies very slowly with angle of 
incidence, even up to fairly large angles. The increase in 
reflectance is only 5% at 30 degrees. This is because rp 


falls as r, rises for angles less than the Brewster angle. 


S 
The results of Cox and Munk were also obtained using 
the same assumption as made here, i.e. that the reflectance 


of the water surface was constant and equal to the normal 


incidence value for all the angles encountered. 
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D. LAMBERTIAN REFLECTION 

In order to compare calculations of predicted returned 
power to actual measured values, it is useful to relate 
glitter reflection to Lambertian reflection, where scatter- 
ing is completely diffuse, as rom a white sheet of paper. 
However the comparison will be made to a Lambertian type 
reflecting surface which reflects the same fraction of the 
total power on it as the reflectance for normal incidence 
from water, i.e. a fraction .0204 of the incident light is 
returned, when integrated over all angles within a hemis- 
phere, 

In general, for Lambertian reflection, the power per 
unit solid angle, I,,;, returned at perpendicular incidence, 
is given by 

Io, = P/n (13) 
where P is the total power reflected from the surface, 
integrated over all angles, 

In the case of an idealized Lambertian water surface, 
the reflected power is given by 

P=rP, 
where, as before, r = .0204 for the water surface. 

The flux at the detector, the power returned per unit 

area, would then be given by 
an eile) ene (14) 
This will be called the "Idealized Lambertian return flux". 


To summarize, it is the flux returning at normal incidence 
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from a Lambertian surface that reflects, integrated over 
all angles, a total fraction of the incident light equal to 


the reflectance of water for normal incidence (.0204). 


E. EFFECTIVE REFLECTANCE 
In order to characterize the optical state of the rough 
water surface, it is useful to define a quantity, pore, the 
effective reflectance of the water surface, as 
Return flux from glitter 


Tie fn a) Seer 
Idealized Lambertian return flux 


of wale a ; (15) 


This quantity can be measured directly in terms of the 
ratio of laser radar return signals from a water surface to 
the signals from a Lambertian surface such as a white paper 
target, often at a different range from that of the water. 
This is diseussed in detail later in this thesis, but to 
avoid confusion it should be mentioned here that the ranges 
of the water and the Lambertian target enter in that case, 
as well as the actual reflectance of the Lambertian surface 
used as a reference, In turn, the anticipated signal from 
a given model of a laser radar altimeter can be estimated, 
using a knowledge of the probable range of pore to be 
expected. 

In terms of the model discussed in section A above, an 
idealized, or predicted, value of Poff would be given by 


equation (15) after substituting expressions (10) and (14), 


Sut 


ere 


p 08 

et ft ~~ ian 
Fou 
[H,] 

or Peff = See e SOa:p) 
26,5, 


This will be called the peeenieen Pests The quantity Bei) 
is the magnitude of the Hermite polynomial for the two- 
dimensional Gram-Charlier distribution at normal incidence. 
Dall had the numerical value of 1.11 in all eases encoun- 
tered here. The quantity pyr-¢ expresses the optical pro- 
perties of a rough water surface, in terms of the standard 
deviations, o, and o,, of the glitter pattern. It is 
particularly useful because it does not depend on range and 
depends only on the optical state of the rough water sur- 
face. 

The functional relationship expressed by equation (16) 
made it seem desirable to plot directly measured values of 
Pafg as a function of G2 CuiGeye The results of such a 


plot are discussed later in this thesis. The measured 


values of pare did prove to vary linearly with 1/(2¢,,¢,). 
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ITI. EXPERIMENTAL PROCEDURE 


A. FIELD SITE 

Experiments to measure laser reflectance as a function 
of glitter intensity profile were conducted from Parrott's 
Ferry Bridge over the New Melones reservoir near Columbia, 
California, following several earlier preliminary experi- 
ments at that site and at the Dumbarton bridge over San 
Franeisco Bay. The Parrott's Ferry Bridge was the most 
favorable of many sites considered as it offered a fairly 
large distance to the water (109 feet at the time of the 
last experiments) and had a pedestrian walkway with no 
obstructions beneath. Heavy automobile traffie and inade- 
quate pedestrian space, as well as obstructions beneath, 
made such sites as the Golden Gate Bridge and other bridges 
in the San Francisco Bay area unsuitable. It is hoped that 
funding for logistie support will be forthcoming so that 
these experiments can be continued from an aerial platform 
over open ocean, 

The data reported here were all collected at the Par- 
rott's Ferry Bridge on the night of 24 November, 1986, a 
clear night with no moon present. The wind speed varied 
from approximately 3 miles per hour to 12 miles per hour, 
Wind direction varied less than 30 degrees. Due to the 
short fetch in the reservoir, varying from 1/4 mile to 1 
mile, depending on the direction, only high frequency waves 


and ripples were observed, It is anticipated that results 
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with these waves will scale up to ocean waves observed from 


higher altitudes. 


B. LASER RADAR ALTIMETER SYSTEM 

The laser system used was a working model of a laser 
radar altimeter, being developed by a companion project. 
The transmitter consisted of a 0.905 micrometer gallium 
arsenide diode laser driven by a pulser having a pulse 
width of 160 nanoseconds at half height. The transmitter 
output beam divergence was 1.4 degrees by 2.8 degrees. A 
silicon avalanche photodiode detector with an 8 mm focal 
length lens was used to receive the reflected laser signal. 
The aperture field of view was circular with a divergence 
angle of 3.6 degrees, totally encompassing the area illumi- 
nated by the transmitter, A 10 nanometer bandwidth multi- 
layer film filter was ineluded in the detector optics to 
eliminate noise due to sunlight, for daylight operation, 
The laser system was mounted on an arm extending 3 feet 
from the bridge railing. There were no bridge supports or 
obstructions near enough to produce unwanted reflections. 
The laser was 109 feet above the water. 

The detector output was continuously displayed on a 
portable oscilloscope. Because the pulse length of the 
returned laser pulses was only 160 nanoseconds, it was 
necessary to stretch the received pulses electronically and 
assemble them to form a continuous-wave analog signal to 


permit recording. By the Nyquist theorem, the maximum 


34 


component frequency of this signal is limited to 500 Hertz 
by the interpulse interval of one millisecond, It was thus 
possible to record this signal on a frequency-modulated 
analog tape recorder, a Hewlett-Packard HP-3960A recorder, 
which had a signal bandwidth of de to 20 kilohertz. The 
signals were recorded simultaneously with recording of the 
video pictures of the glitter pattern. Synchronization of 
the video recordings and laser return signal recordings was 
accomplished by recording the same voice track on both 
recorders. 

The laser return signal magnitude was later data-pro- 
cessed in the laboratory with a DATA Precision Corp. DATA- 
6000 waveform analyzer to yield a mean laser return signal. 
As the laser return signals were played back from the 
recorder into the analyzer, 16 sequential waveforms, each 
consisting of 512 data points taken at 1 msec intervals, 
were averaged. This process was repeated 5 times for each 
run and these were averaged to obtain a mean value. Suc- 
cessive averagings of slightly different regions within the 
ten seconds over which the video signals were averaged 
yielded a standard deviation for the magnitude of the mean 


Signal. 


C. LASER REFERENCE MEASUREMENTS 
An effective Lambertian water reflectance signal was 
obtained by means of measurements of the laser return 


signal from a white sheet of paper. These measurements 
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were made in the field shortly after the measurements of 
water reflectance signal. The white paper sheet was lo- 
cated at a distance of 39 feet. At this distance the 
transmitted laser pattern was still fully within the boun- 
daries of the sheet of paper. The same paper was separate- 
ly calibrated in the laboratory and found to reflect 0.73 
times as much light at normal incidence as a perfect Lam- 
bertian surface, This value is representative of most 
white paper surfaces, The calibration procedure is dis- 
cussed in Appendix B. The effective Lambertian water re- 
turn signal was then obtained by multiplying the signal 
from the white sheet of paper at 39 feet by the ratios 
(.0204/0.73) x (39/109)? . 

The factor (.0204/0.73) takes account of the difference in 
reflectance of the white paper sheet and the idealized 
Lambertian water surface. The factor (39/109)% takes ac- 
count of the different distances of the white sheet and the 
water, 39 and 109 feet respectively. The inverse square 
law applies to these distances. 

The effective reflectance coefficient, pore, as defined 
in section II E, was then calculated by dividing the mean 
laser water return signal by the effective Lambertian water 
return signal. i.e 

Mean laser return signal 


Poff SO a ea aaa e 
Idealized Lambertian water return signal 


Se 


D. GLITTER PATTERN MEASUREMENT 

Glitter patterns were produced by illuminating the 
water surface with a quartz-halogen lamp of a type used as 
a light source in overhead classroom projectors. This pro- 
duced a nearly uniformly illuminated patch on the water 
with a total beam angular divergence of about 60 degrees, 
considerably wider than any of the observed glitter pat- 
terns. The uniformity of the illumination was verified in 
the laboratory in advance. The lamp was rated at 360 watts 
at 54 volts. Power for this lamp was provided by a small 
rotary converter generating nominal 120 volts, 60 Hertz 
power up to 500 watts, and driven by 12 volts de from an 
automobile-type storage battery. The lamp voltage was var- 
ied, as needed, by means of a variac. The video camera 
was an RCA type TC20595C with a Vidicon tube. It was oper- 
ated on 12 volts de directly from a second automobile 
storage battery. The camera performance was very stable as 
it generated its own line-scan frequency with a quartz 
oscillator. The HP3960A FM data recorder for the laser 
radar signals also operated with 120 volt 60 Hertz power, 
obtained from this same battery through a frequency con- 
trolled square wave chopper converter. The video images 
from the RCA camera were recorded on a portable 8mm video 
tape recorder that was part of an Olympus VX-801-KU 8mm 
video camera system. This system also displayed a contin- 


uous monitor image of the scene as viewed by the RCA 
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camera. This allowed the brightness level of the illumina- 
tion to be adjusted for proper recording and also permitted 
alignment of the camera with the major or minor axis of the 
elliptic glitter pattern. The 8 mm tape system was very 
compact and operated on its own internal batteries. As 
measured in the laboratory, the resolution of this combined 
system exceeded that of the usual home VHS format by about 
20%. 

The lamp and video camera were mounted on an arm exten- 
ding from the bridge railing at a point about 6 feet from 
the laser altimeter. The video camera was oriented first 
along, then across, the wind direction, so that sequential- 
ly the major and then the minor axis of the elliptical 
glitter pattern would be aligned horizontally in the re- 
corded video image. The 16 mm camera lens was focused at 
infinity with an aperture setting of f/1.6 . The 16 mm 
focal length gave a wide enough field of view to include 
the complete glitter pattern in every case, Recording 
sequences of 15 to 30 seconds were made at several differ- 
ent illumination levels in order to be sure to avoid satur- 
ating the camera, 

Later, in the laboratory, the 8 mm video tape recor- 
dings were copied into a Sony Superbeta VCR which had more 
flexible reproducing features than the 8 mm field unit, 
The resolution of this VCR was previously measured and 
found to exceed that of the usual home VHS format by 50%. 


No detectable degradation occurred in this copying process. 
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The output of the video recorder was sent to a Tek- 
tronix 468 digital oscilloscope. Successive single hori- 
zontal TV scan lines, lying along, or perpendicular to, the 
glitter pattern major axis, were displayed on the scope 
using the "B sweep delayed" mode. The waveform displayed 
was a time-varying signal proportional to the intensity 
distribution along the major or minor glitter pattern el- 
liptie axis. The oscilloscope was then set up to digitize 
and average 256 such sweeps and to store the accumulated 
waveform, The real time period represented by the 256 
sweeps was 256/30 or 8.5 seconds. The computer program in 
appendix A was written to allow the Tektronix 468 to inter- 
face with an IBM PC/XT through an IEEE488 interface bus 
bref. 7]. With the computer as controller for the opera- 
tion, the oscilloscope digitized and transmitted the waves 
form data to the computer, The digital information was 
then stored in memory or on disk for further processing. 

The Tektronix oscilloscope was then adjusted so that it 
displayed a TV line through an image of the prevailing 
water surface outside the glitter pattern. Sweep averaging 
and transmission of a waveform to the computer memory were 
again carried out. This provided a background signal which 
was subtracted from the glitter waveform in the computer. 
The computer was then programmed to plot the glitter pat- 
tern profile on an Epson F X80 dot matrix printer. Figures 


3 and 4 show the resultant intensity profiles along the 
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major and minor elliptic axes from a representative run. 
Gram-Charlier standard deviation values o, and o, were then 
measured graphically from these profiles by measuring the 
half width at a height equal to 0.659/1.11 times the peak 


value, as discussed earlier, 


42 


IV. RESULTS 


The measured values of pire, the ratio of observed 
laser radar return to that expected for Lambertian reflec- 
tion, were plotted as a function of the corresponding 
Pabues ofgW/(2Zc76,),9in Figure 5, where ¢o 


and o. were the 


u Cc 


measured upwind and crosswind Gram-Charlier standard devia- 
tion values for the glitter patterns. The data values are 
tabulated in Appendix C. Use of the function 1/(2¢,o,) for 
the abscissa was suggested by the model discussed in sec- 
tion II. This functional form organized the experimental 
points so as to fall fairly well along a straight line. 
The solid straight line was a least squares solution for 
the data points taken here, shown as the circles in that 
figure, The correlation coefficient was r = 0.96 for the 
points relative to this line. This line had a slope of 
0.92 and an intercept of 44 on the ordinate axis. 

This empirical relationship provided a basis on which 
to predict the expected reflectance of a pencil beam from a 
rough water surface on the basis of only the remote sensing 
measurement of the glitter pattern profile. This made 
possible the field evaluation of working models of laser 
radar altimeters with the help of a simple technique. The 
functional relationship, and also the body of data were the 
first of this sort in this field. Without this technique, 
the field performance of a given laser altimeter system 


would require expensive concurrent logistic support to 
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Figure 5. Observed effective reflectance relative 
to Lambertian, pore, as a function of 1/(2¢,,¢,), 
plotted as circles, The solid line is a east 
squares fit of a straight line to these points. Two 
data points plotted as crosses are based on Petri's 
data [Ref. 2]. The dashed line represents an ideal 
model proposed in this thesis. 
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determine the roughness state of the water. The latter was 
likely to be inferred from meteorological and oceanographic 
data, and to be only the deduced probable steady state 
condition, whereas the technique reported here directly 
measured the momentary optical state of the ocean surface, 

The numerical values of pure were of immediate utility. 
The lowest value, represented by the intercept on the y 
axis, was 44, Although extrapolation to the axis had con- 
siderable uncertainty, it did imply that for the worst case 
the ratio of reflectance to that for Lambertian was con- 
siderably above unity, at least by somewhat over one order 
of magnitude. Without evidence such as this, a working 
altimeter would presumably have had to be designed to cope 
with a case where a unity ratio to Lambertian might have 
occurred. This factor of 44, or somewhat over one order. of 
magnitude increase in signal, represented a large reduction 
in the design requirement to be sure of a return signal. 
This had immediate importance in the ability to achieve the 
design requirements. 

That the observed signals were directly relatable to 
reflection from a Lambertian surface made it possible to do 
most of the testing work in the laboratory. Various evolu- 
tions of laser radars could be tested with a Lambertian 
surface such as a white sheet of paper. The signals could 
then be scaled to those for a water surface at any required 
distance. The range of field signal return values could 


then be evaluated in terms of Figure 5. 
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The model proposed in Section IIB, which predicted that 
Poff = LH, l/c some required a straight line through the 
origin with slope [H,] = I.1! in Figure 5.” This=is Tepmes 
sented by the dashed line in that Figure. An alternative 
prediction is for the slope to be [H,J]/4 , depending on the 
view taken regarding the solid angle of the return scatter- 
ing cone, as described in Appendix D. Although neither 
alternative represented a good fit, the data points did lie 
along a straight line, with slope not far from 1.11. The 
existence of a theoretical model that agreed with the 
observed behavior in more detail would have been desirable, 
but it was not necessary. The present experiment had 
achieved the stated objective of providing a means of 
determining the optical state of the rough water by a 
remote sensing technique, where none existed before this 
work, The result was semiempirical, but was useful. 

The fact that the proposed model did not predict the 
observed behavior could be viewed as an interesting chal- 
lenge. That the observed effective reflectance was larger 
than expected, by somewhat less than a factor of four, and 
that the straight line had an intercept, indicated that 
this was a problem unlikely to be related to miscellaneous 
small increases or losses of light. Some of the immediate, 
ad hoc, explanations of this discrepancy that come to mind, 
seemed to predict less reflectance, rather than more, For 
example, if the variation of specular reflectance with 


angle were invoked for light reflected at large angle of 
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incidence into the far wings of the Gram-Charlier distribu- 
tion, then less light should have returned along the axis, 
rather than more, 

Another effect producing deviation toward too little 
reflection was that the laser "pencil" beam had a finite 
divergence of 1.4 by 2.8 degrees. Thus the reflectance was 
really the average over the curved top of the Gram-Charlier 
distribution. Thus the average observed would have been 
slightly less than the peak value. Correction for this is 
planned in subsequent work, but this correction will pro- 
duce an increase in the peak reflectance, rather than a 
decrease, as needed to fit the model. 

AMPOsSSIOllLIty. In the right direction, in that it would 
have tended toward producing the higher than expected val- 
ues of reflectance, was that some spray might have begun to 
develop as the wind increased. This would have correspond- 
Paeeomne repionwor "small values of 91/(20,90,). The drop- 
lets would have acted as retroreflectors and increased the 
returned signal along, or near, the laser beam axis. This 
effect, called the "glory", is commonly seen in looking 
down from an aircraft in the direction of the shadow of the 
aircraft. However, this explanation seemed quite unlikely 
in view of the low sea state encountered on the reservoir. 

As to the major factor of four, it appeared that there 
was no real reason to omit it. The belief that it possibly 
lay in misinterpretation of the solid angle of the return 


cone of light had its only real motivation in that the 
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alternative choice led to a factor of four. It seemed more 
likely that the trouble was fundamental, perhaps lying in 
the assumptions made by Cox and Munk, and others, that the 
light was reflected by infinitesimal flat facets, rather 
than by a statistical average over focusing curved sur- 
faces, where the surfaces fortuitously curved to be retro- 
reflectors contributed disproportionally high return sig- 
nal, Such an effect seemed likely to be much more pro- 
nounced with a nearly point source at short range, as used 
here, than for the case of the sun as an extended source at 
infinite range, as used by Cox and Munk. This point is 
expected to receive more attention in the succeeding work 


on this project. 
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V. COMPARISON WITH OTHER RESULTS 


Direct comparison of the results reported here with 
those of others was possible only for one bit of frag- 
mentary data. No other measurements of laser reflectance as 
a function of glitter pattern width were Known to the 
author, except for a single diagram in the article by Petri 
[Ref. 2]. The values of glitter pattern width were not 
quoted there, as the article was primarily concerned with 
measuring the peak reflectance as a function of wind velo- 
city. However, a few sample glitter pattern profiles were 
shown in one figure. That figure is reproduced here as 
Figure 6. It was possible to measure the sigma values 
graphically from these curves, although the curve shapes 
were poor, Those curves were obtained by a laser scanning 
in a vertical plane, with the plane rotating slowly through 
360 degrees about a vertical axis. The scans were averaged 
so that no information was retained on the relative magni- 
tude of upwind and crosswind sigma values. Two of these 
curves yielded sigma values that could be related to the 
results in this thesis. The others were too poorly re- 
solved or they corresponded to reflectivity above or below 
the range measured in this work. Two pairs of these curves 
were for essentially the same windspeed so that an estimate 
of the internal consistency could be made. These pairs 
differed by 12% and 14% for the peak value. Points for the 


two utilizable widths were plotted with the symbol + on 


49 


NOTE: Each normalized curve 


0.7 represents the average of all . 
18 azimuth positions 27 Jun - 3 mph (estimated) 


"Exceptions (5 minute averages 
in one azimuth position) 


All data was collected with a 
0.6 scanner speed of 16 scans/sec 


and a field of view of 24 
milliradians 


0.5 
28 Jun - 7.4 mph 


0.4 
*29 Jun - 9.6 mph 
0.3 
4 Oct - 9.6 mph 
*29 Jun - 10.4 moh 
29 Jun - 10.8 mon 
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Figure 6. Normalized average return signal versus 
laser/receiver angle as reported by Petri [Ref. 2]. 
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Figure 5 of this thesis. They fell close to the data points 
obtained in this work, 

Although the results obtained in the work reported here 
were intended to provide a means of evaluating laser radar 
altimeter performance without the.necessity of knowing the 
wind speed, it was interesting to relate these results to 
data taken by others, where the wind speed was measured, 
This provided some confirmation of the validity of the 
results reported in this thesis. 

An indirect comparison, where an inferred wind speed 
allowed intercomparison with related experimental data in 
the literature, was possible by combining the results of 
Petri [Ref. 2} with the results of Cox and Munk [Ref. 3]. 
Petri reported laser reflectance data in terms of wind 
speed for waves under the Chesapeake Bay bridge near Anna- 
polis, Md. Figure 7 reproduces the results of Petri for 
laser return as a function of wind speed, with his results 
plotted as the crosses. No theoretical model was offered 
in that article, 

Cox and Munk gave a relationship between their measured 
Wind speed and the mean square upwind and crosswind sigma 
values for the glitter patterns, where W was the wind speed 
in meters per second, as 

moe GOS toed x 10e° We.004 . 

W has been called the effective wind speed here, as it was 
only a means to cross-connect to the data of Petri. Using 


the above expression, effective wind speeds were calculated 
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Figure 7. Plot of observed effective reflectance rela- 
tive to Lambertian, p.¢g¢, versus wind speed. Cireles 
represent reflectance da a obtained in this work plot- 
ted versus inferred wind speed. Crosses represent ref- 
lectance data collected by Petri plotted versus mea- 
sured wind speed [Ref. 2]. 


DZ 


for our data. The corresponding laser returns are plotted 
as the circles in Figure 7. It can be seen that the data 
points from both sources cluster into a broad band in the 
same general region. 

The approximate agreement of the data reported here 
with the combined results of Petri and Cox and Munk gave 
some credence to the ability to relate laser radar return 
to the width of glitter patterns. 

It should be pointed out again that the purpose of this 
work was to obtain a measure of expected laser radar return 
signal from rough water by remote sensing techniques and 
without a knowledge of wind speed. The effective wind 
speed deduced from the expression of Cox and Munk was 
useful only to permit the cross-connection to Petri's data. 
These wind speeds should be considered here only as a means 
of obtaining order of magnitude confirmation of the optical 
results. Intercomparison was somewhat doubtful because the 
Cox and Munk data were for open ocean waves, with wind 
Speeds measured at 9 and 41 feet above the water. The 
Petri results were for waves under the Chesapeake Bay 
bridge near Annapolis, Md, with wind speed measured at 60 
feet above the water. In spite of these uncertainties in 
wind velocity it was interesting that utilization of these 
results yielded approximate agreement between the data re- 


ported here and that of Petri. 
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VI. CONCLUSION 


This work showed the feasibility of a remote sensing 
technique for determining the expected magnitude of laser 
radar return from a rough water surface by measurement of 
the size of a simultaneously measured elliptical glitter 
pattern. The technique did not require any additional 
Knowledge of the water surface beyond optical measurements 
made at the location of the laser radar altimeter. No 
previous such technique had been in existence. The effec- 
tive laser reflectance, pyre, proved to vary approximately 
linearly with 1/(2¢,¢,), where o, and ¢, are the Gramm 
Charlier standard deviations for the upwind and crosswind 
glitter pattern intensity profiles, respectively. pur, is 
the ratio of the laser radar return signal to the ‘laser 
radar return signal to be expected from the water if it 
were to act as a Lambertian, diffuse reflecting, surface 
but with an overall reflectance coefficient equal to that 
for normal ineidence for water. A least squares fit to a 
straight line gave a slope of 0.92 and a y intercept of 44, 
with a correlation factor, r = 0.96 . 

Plotting pure as a function of 1/(26,¢,) was suggested 
by the model proposed in this thesis. That model indicated 
that a straight line for pore as a function of 1/(2¢) ae 
should pass through the origin with a slope of [H,]/4, or 
perhaps a slope of Distal where [H,] = 1.11 , under a dif- 


ferent, controversial, interpretation. The solution to the 
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discrepancy between the observed data and the suggested 
model is a problem that remains for further work. 

That the straight line did not go through the origin, 
but instead had an intercept of 44, provided important 
information. This said that the worst case of very rough 
water should provide a return signal 44 times the magnitude 
of that from a water Lambertian surface, Although such 
extrapolation completely to zero is risky, it implied that 
for very rough water the smallest laser radar return is 
somewhat more than one order of magnitude larger than would 
have had to be assumed ie the absence of this information. 
This would permit the design of a laser radar altimeter 
with an order of magnitude less power than would have to be 
assumed otherwise. 

Finally, the existence of these results provides a firm 
relationship between the reflectance of a rough sea surface 
and that of a Lambertian surface, such as a large white 
surface at limited range in the laboratory. This reduces 
the need for field tests of various working models of laser 


radar altimeters, 


DD 


APPENDIX A 


COMPUTER PROGRAM 


DIM Y(520) ,DTA(520) , BACGND(520 
DEF SEG = &HDOOO 


KE ere : 
INIT = 0 "ASSIGN OFFSET ADDRESS OF Via 
ADDR% = 21 “CONTROESER ADDREoS 


LEVEL% = 0 

CALL INIT(ADDR%, LEVEL?) 

TRANSMIT= 3 *KECEIVE = 6 “"OFFSED ADDR EcorS 
CMDS ="MLA TALK 3" 'MY LISTEN ADDRESS 


'DEVICES' TALK ADDRESS 
CALL TRANSMIT (CMD$,STATUSS) 


PRINT "TRANSMIT STATUS =" "7olATUss 


R$ = SPACES(255) ‘SET LENGTH OF RS>= LENGTH 
R1S = SPACE$S(255): R2$ = SPACES(255) 
PRINT "Set up for data - push any key when ready " 


VS=INKEYS:iF VS="" THEN 150 

CLS 'OF THE INCOMING DATA 
CALL RECEIVE(R$, LENGTHS, STATUS$) 

PRINT "RECEIVE STATUS = ",STATUS3 

CALL RECEIVE (R1$,LENGTH2%,STAT2%) 

PRINT "RECEIVE STAT2 = ",STAT2% 


CALL RECEIVE(R2$,LENGTH3%,STAT3%) 
PRINT "RECEIVE STAT3 ",STAT3% 

IF STATUS%<>0 THEN PRINT'RECEIVE ERROR TRY AGAIN" ;GOTO 
100 

FOR I = 178 TO 255 

DTA (I-177) =ASC (MID$(r$, i) ) 

NEXT I'270 PRINT R$,R1$,R2$ 

FOR I=1 TO 255 

DTA(I+78) = ASC(MID$(R1$,i00 
NEXT I 

FOR I= 1 TO 180 

DTA (333) =ASC(MID$(R2$,1)) 

NEXT I 

FOR I=1 TO 520 STEP 10 

FOR J = 1 TO 10 


PRINT USING " ##H2" > DTA (I+J—1) 7 

REM LPRINT USING" HEU sdta (I+J— L) 3 ' 380, NEXr J 
PRINT 

REM LIST 

NEXT I 

INPUT "DO YOU WANT THE DATA RECORDED? Y/N", ANSS$ 
IF ANS$ = "n" OR ANSS$ = "N", THEN GOTO 580 


INPUT “INPUT FILE NAME", FILENS 
IPUT "input header",A$ 
OPEN FILENS$ FOR OUTPUT AS #1 
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470 PRINT #1,A$ 

480 PRINT #1, MIDS(R$,1,80) 

490 PRINT #1, MIDS$(R$,81,80) 

BOG PRINT #1, MIDS(RS, 161,13) 

510 PRINT #1, DATES, TIMES 

520 FOR I=1 TO 520 STEP 10 

530 FOR J=1 TO 10 

540 PRINT #1, USING " ###" > DTA (I+J-1) ; 
550 NEXT J 

560 PRINT #1," " 


570 NEXT I 

580 R$ = SPACE$(255) 'SET LENGTH OF RS$>=LENGTH 
590 R1$ = SPACE$(255): R2$ = SPACES(255) 

600 'OF THE INCOMING DATA 


610 PRINT “Set up for background - push any key when ready " 
620 VS=INKEYS: IF VS=""" THEN 620 


oe. CLS 

640 CALL TRANSMIT (CMDS,STATUS%) 

650 PRINT "TRANSMIT STATUS = ",STATUS% 
Seu, CALL RECEIVE(RS, LENGTHS, STATUS?) 
670 PRINT “RECEIVE STATUS = ", STATUS% 
e50) CALL RECEIVE (R15, LENGTH2%,STAT23) 
e20 PRINT “RECEIVE STAT2 = ",STAT2% 
700 CALL RECEIVE (R2$,LENGTH3%,STAT33%) 
fo PRINT “RECEIVE STAT3 = ",STAT3% 


720 IF STATUS%<>0 THEN PRINT"ERROR READING DATA ":GOTO 580 
722 IF STAT2%<>0 THEN PRINT"ERROR":GOTO 580 
724 IF STAT3%<>0 THEN PRINT"ERROR":GOTO 580 
730 FOR I = 178 TO 255 

740 BACGND(I-177) =ASC(MID$(R$,I)) 

750 NEXT I 

760 PRINT R$,R1$,R2$ 

770 FOR I=1 TO 255 

780 BACGND(I+78) = ASC(MID$(R1$,1I)) 

790 NEXT I 

BeGmEOR I = 1 TO 180 

810 BACGND (I+333)=ASC(MID$(R2$,I)) 

820 NEXT I 

830 FOR I=1 TO 520 STEP 10 

B40 FOR J = TO 10 


850 PRINT USING " ###" > BACGND(I+J-1) ; 

860 REM LPRINT USING " ###" > bacgnd(I+J-1) ; 

870 NEXT J 

880 PRINT 

890 REM LIST 

900 NEXT I 

910 INPUT "DO YOU WANT THE DATA RECORDED? Y/N", ANSS 
920 IF ANSS = 'n' OR ANSS = "N", THEN GOTO 1080 


930 INPUT "INPUT FILE NAME", FILENS 
940 INPUT "input header",AS 

950 CLOSE #1 

960 OPEN FILENS FOR OUTPUT AS #1 
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970 PRINT #1,AS$ 

980 PRINT #1, MIDS(RS$,1,80) 

990 PRINT #1, MIDS$(RS$,81,80) 

1000 PRINT #1, MIDS( RS, ter, Tey 

1010 PRENT #71, DATES, TEMES 

1020 FOR =f TO 520 Siler ere 

1030 FOR J=1 TO 10 

1040 PRINT #1, USING " ###" > BACGND(I+J-1) ; 

1050 ° NEXT J 

1O0GO2ERINT #17) 

1070 NEXT I'1075 YMAX=-512 

1076 YMIN=512 

1080 FOR I=1 TO 520 

1090 Y¥(I)=DTA(I) - BACGND(I) 

1091 EF Y¥ (21) <YMEINSEHEN YMIN=yY ia 

1092 IF Y(I)>YMAX THEN YMAX=Y (1) 

TLOOmiE XT sl 

1120 FOR IJ=1 TO 520 STEP 10 

1130 FOR J=12 foro 

1140 PRINT USING " HHHHW SY (I+J-1) ; 

1150 WExX Toe! 

1160 PRINT 

117 Oe NE 

1180 XMIN=100 

1190 XMAX=364 

1500 YWMIN=YMIN-.2* (YMAX-YMIN) 

1510 YWMAX=YMAX+.15* (YMAX-YMIN) 

1520 XWMIN=XMIN-.2* (XMAX-XMIN) 

1530 XWMAX=XMAX+.15* (XMAX-XMIN) 

1540 SCREEN 2,1 

15507 few Ore) (Geo 1 oo) 

1560 WINDOW (XWMIN, YWMIN) - (XWMAX, YWMAX) 

1580 YLABS="RELATIVE INTENSITY" 

1590 XLABS=""POSITION IN GLITTER PATTERN" 

L6OOREES 

1620 IXX=PMAP (XMAX, 0) 

1630 IXM=PMAP (XMIN,0) 

1640 IYX=PMAP(YMAX,1) 

1650 IYM=PMAP(YMIN,1) 

1670 DRAW"BM=IXM; ,=IYM;'"': DRAW"M=IXM; ,=IYX; :DRAW"M=IXX;, 
=Ly x" 

1680 DRAW"M=IXX; ,=IYM;": DRAW"M=IXM; ,=I1YM;" 

1681 EOCATE 235525 

1682 PRINT XLABS 

1700 FOR X=XMIN TO XMAX 

1710 IX=PMAP(x,0) 

1720 IY=PMAP(y(X),1) 

1730 IF X=XMIN THEN DRAW"'BM=I1X; ,=1Y;" 

1740 DRAW"M=I1X; ,=1Y;" 

17 SO Shiner +x 

1751 LOCATE 6,4: PRINZ PORELATLIVE” 

1790°WIDTH ipe eo. 
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1800 
2000 
2005 
Z0 20 
3000 


LPRINT CHRS$ (27) "A"CHRS (7) 
VS=INKEYS:IF VS="" THEN 2000 
LPRINT CHR$(27)"A"CHR$ (12) 
SCREEN 0 

KEY ON 
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APPENDIX B 


LAMBERTIAN REFLECTANCE CALIBRATION 

Calibration of the Lambertian reflectance coefficient 
for the white paper used as a reference in the field was 
obtained by measuring the ratio of laser radar return 
signals from a small segment of area, Ay, of the white 
paper target, to the return signals for specular reflection 
from a plate glass window, 

The Lambertian reflectance, r,;, of the the white paper 
is given by 
Tale Py ae 


Ane Tea The 


fy = 


where 


rs = Reflectance of the plate glass window, for two 
surfaces = .0826 for n = 1.51, 


P, = Laser signal from a small area Lambertian 
target, 


P; = Laser signal from a plate glass window, 

R; = Distance to the Lambertian target, 

R, = Distance to the plate glass window, and 

A, = Area of the white Lambertian target, small 
enough to lie within the central flat 
illumination region of the laser transmitted 
beam, 

The calibration described above was independent of the 


Lambertian properties of the calibrating target of white 


paper. That is, it did not matter if the target did not 
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follow the Lambertian behavior of brightness being propor- 
tional to the cosine of the angle from the normal. This 
was because the calibration technique involved only mea- 
surements made with the laser radar and hence with a small 
acceptance angle near normal incidence, 

The flatness of the plate glass window used in the 
calibration was measured by determining the position of the 
reflected image of an incandescent lamp source. This per- 
mitted calculating the radius of curvature, The glass was 
then treated as a curved mirror, although this correction 
was very small. 

The value of Lambertian reflectance obtained was 0.73 . 
This is comparable to values of about 0.75 given in many 
tables of values for the Lambertian reflectance of white 


paper. 
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APPENDIX C 


TABLE OF DATA 


16 


17 


Glitter Profile Effective Laser Inferred 

Standard Deviations Reflectance Wind Speed 
(Pere) 

Upwind Crosswind Predicted Observed 

(o,,) (a,) 
~340 .094 17.4 46.8 a 4 
.224 ~121 20.9 47.6 7.18 
2175 ~141 22.9 42.0 2.60 
~134 .058 71.4 90.9 2.41 
A422 026 175 150 Wie: 
.163 .024 142 182 2.90 
.081 030 228 213 L207 
.084 .034 194 217 0:92 
-2io 064 390 89.8 9.69 
181 058 92.9 Dead 4.04 
~241 120 EE) 7 73.0 8.45 
135 118 34.8 88.3 4.80 
~167 .110 30.2 P26 3 4.43 
~150 .100 37.0 86.9 3.58 
~126 046 95.8 116 2.39 
~217 .201 Ar 43.5 10.4 

301 096 192 42.3 11.0 
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APPENDIX D 


LIGHT SCATTERING CONE 

In calculating the expected laser reflectance from a 
rough water surface, it is necessary to evaluate the size 
of the cone into which the surface reflects an incident 
narrow laser beam, However, from geometrical opties, this 
cone does not have the same included solid angle as the 
cone measured by observing the glitter pattern obtained 
with an imager located near a point source illuminator, In 
the latter case, the light reflects at exactly normal 
incidence from the facets on the water surface, and the 
outer limits of the glitter pattern (or sigma points, since 
a definite outer limit does not strictly exist) occur at a 
definite angle of inclination of the water surface, How- 
ever for a vertical pencil beam, the reflected light re- 
turns into a cone that has twice the angular spread in both 
x and y directions, and hence a solid angle that is four 
times that of the glitter pattern observed at the source 
with a point source illuminator. This is because the 
surface which the pencil beam encounters contains facets 
that are inclined at all the angles which the glitter 
pattern indicates are in existence, if averaged over a 
significant period of time. Thus, for example, some of the 


pencil beam light will encounter surface facets that are 
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inclined at the angle corresponding to the sigma point of 
the glitter cone. Because this light is now not at normal 
incidence to the facet it is encountering, the angle of 
reflection will be equal to the angle of incidence and the 
total angle of deviation will be twice the angle of inci- 
dence, Since the angle, at the glitter pattern position 
corresponding to the sigma point, represents the angle of 
inclination of a facet surface, the angle of deviation from 
the vertical of the light from a vertically downward pencil 
beam will be at an angle equal to twice the angle for the 
glitter pattern sigma point. Since the deflection can be 
in either the x or y direction, the solid angle of the cone 
into which the light from the pencil beam is reflected 
would be expected to have four times the solid angle of the 
glitter pattern, 

Unfortunately, the results obtained with the above 
interpretation lead to a predicted laser return that is too 
small by approximately a factor of four, when compared to 
the observed laser return signal. At this writing, the 
explanation of this anomaly is not at hand. In order to 
lead to a reasonable fit to the observed laser return, it 
has been necessary to assume that the cone into which the 
reflected light radiates is equal in size to the glitter 


pattern cone. 
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